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Abstract. The Enskog-Landau kinetic equation is considered to describe non-equilibrium processes of a 
mixture of charged hard spheres. This equation has been obtained in our previous papers by means of the 
non-equilibrium statistical operator method. The normal solution of this kinetic equation found in the first 
approximation using the standard Chapman-Enskog method is given. On the basis of the found solution 
the flows and transport coefficients have been calculated. All transport coefficients for multicomponent 
mixture of spherical Coulomb particles are presented analytically for the first time. Numerical calculations 
of thermal conductivity and thermal diffusion coefficient are performed for some specific mixtures of noble 
gases of high density. We compare the calculations with those ones for point-like neutral and charged 
particles. 
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Construction of kinetic equations for dense gases and 
plasma is one of the most important problems in the ki- 
netic theory of classical systems. A sequential kinetic the- 
ory of dense systems does not exist yet. The Enskog- 
Landau kinetic equation has been obtained recently in pj 
to describe transport processes in non-equilibrium system 
of charged hard spheres. This equation has its name due to 
the structure of total collision integral. This integral con- 
tains terms of revised Enskog theory, kinetic mean field 
theory and Landau- like collision integral [0. The influ- 
ence of the last term on system behaviours has been of 
our main interest. This term is caused by the long-range 
interactions in a system. In particular, it was shown P[ 
that in the case of small densities and weak interactions 
the force autocorrelation function and the entire last term 
in the total collision integral converts to the usual Landau 
collision integral for a rarefied plasma. 

The new kinetic equation and its collision integral are 
adequate only for systems which can be modelled by char- 
ged hard spheres. The great credit for this result is shared 
by a choice of interparticle interaction potential in additi- 
ve-like form: hard spheres interaction plus certain "smo- 
oth" long-range part (Coulomb interaction in our case). 
It allowed to avoid divergency at short distances. Unfor- 
tunately, a logarithmic divergency at long distances still 
remains, and to eliminate it one should introduce a cut- 
off radius of integration (like a Debye one). But unlike 
the classical Debye formula, in this case we used a modi- 
fied one, which takes into account particle sizes a. For the 



Enskog-Landau kinetic equation |l| a normal solution has 
been found by means of the standard Chapman-Enskog 
method. A stress tensor 11 and heat flow vector q have 
been obtained as well. Expressions for transport coeffi- 
cients like bulk k and shear rj viscosities and thermal con- 
ductivity A have been derived from structures of l~l and q. 
Numerical calculation of transport coefficients for neutral 
and once-ionized argon shows a good agreement between 
the theory and experimental data. In |||4|] these results 
were generalized to non-stationary non-equilibrium pro- 
cess. Whereas to find the normal solution the Chapman- 
Enskog method || is used in the much morepower- 
ful method of boundary conditions M is used in In 
the limiting case of a stationary non-equilibrium process, 
the results of j| completely convert to those of Q. For 
hydrodynamic description of fast processes it is better to 
use the method of boundary conditions Q . 

Application of the theory to a multicomponent system 
was performed step-by-step. The Enskog-Landau kinetic 
equation for M-component (M ^ 2) mixture of charged 
hard spheres has been proposed in |7]]. Just the same, the 
normal solution, flows and transport coefficients have been 
found by means of the standard Chapman-Enskog method 
for a two-component system only. New transport coeffi- 
cients which appear in multicomponent systems are mu- 
tual diffusion D a @ and thermal diffusion D^, coefficients 
(here a, (3 are mixture indices). Numerical calculation of 
the obtained transport coefficients showed a good agree- 
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merit between the developed theory, experimental data, 
results of other theories and MD simulations. 

In view of identifying of a normal solution of the kinetic 
equation for a multicomponent system of charged parti- 
cles results of |?| are not the first. Namely, it is worthy to 
note that the normal solution of the kinetic equation for 
completely ionized plasma is found in || using the stan- 
dard Chapman-Enskog method in the 3rd approximation. 
But dense systems of finite size particles are consistently 
considered for the first time in our papers []l],||,||,[|,f7]]. 

In this letter we present our solution of the Enskog- 
Landau kinetic equation for a multicomponent mixture of 
charged hard spheres. Similarly to 0, we use the stan- 
dard Chapman-Enskog method. It is known [|| that the 
correction expression for one-particle distribution function 
in the first approximation can be chosen in two different 
ways. According to the first one described in II, the cor- 
rection for one-particle distribution function of cv-kind is 
proportional to d Q - a diffusion thermodynamic force of 
cv-kind only. The second way is proposed in Q. In this 
case the correction is proportional to a certain superposi- 
tion of diffusion thermodynamic forces of all components 
of a mixture. It was shown jl0| that the second method 
gives much better results because, unlike the first method, 
after crossing to linear thermodynamics equations it is 
important that Onsager's reciprocal relations obey. Nev- 
ertheless, in [ij for a two-component system we used the 
first method because the complication (or generalization) 
like in JsJ is essential for two-component systems in some 
cases only. Namely, when the density of particle number 
of some mixture component (or components) is not con- 
served. Such a situation can be realized in gas mixtures 
where chemical reactions between components may take 
place, or in multicomponent mixtures, where transitions 
between states with different internal degree of freedoms 
are possible. 

Let us consider the Enskog-Landau kinetic equation 
for multicomponent mixture of charged hard spheres: 
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here /i denotes the one-particle distribution function, x = 
{r; p} is a set of phase coordinates in a phase space while 
r and p denote the Cartesian coordinate and particle mo- 
mentum, respectively. Collision integral of this equation 
has an additive structure 
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and each term in (|J) is defined like in j?J : the first and sec- 
ond ones are from so-called hard sphere part of interpar- 
ticle interaction - collision integral of the revised Enskog 
theory RET: 
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the third one is caused by taking into account of long- 
range interparticle interaction in the mean field approxi- 
mation KMFT (this term is of the first order in interac- 
tion): 
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and, finally, the third one is the so-called Landau-like col- 
lision term (it is of the second order in interaction): 
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In expressions for 1^ , Ijg , Imf and Jl we used the follow- 
ing designations: b - impact parameter, (3 - an analogue of 
local inversed temperature, e - azimuthal angle of scatter- 
ing, g, a @ - relative velocity of a- and /3-kind particles, g^ 
- two-particle correlation function, m* - reduced mass, 
m a - partial masses of particles, n - total density of par- 
ticles number, n a - partial densities of particles numbers, 
r^2 ~ unit vector along r"^ direction, v' - velocities of 
hard spheres after collision, 9{. . .) - Heaviside unit step 
function. 

In the limit of a system of point-like charged parti- 
cles of low density, the collision integral (13) transforms 
to the usual Landau collision integral Following 
a concept of consistent description of kinetics and hydro- 
dynamics of non-equilibrium processes |fl3, 14, 15|,[l6|,^7|, 
kinetic equation p) should be solved together with local 
conservation laws [H for additive invariants. These addi- 
tive invariants in collisions (or scattering) of charged hard 
spheres are mass (or total density), momentum and total 
energy ||||,[ll],[l2|]. It should be noted that for rarefied sys- 
tems it was sufficient to consider kinetic energy only, while 
in dense systems potential interaction energy is essential 
and cannot be neglected. 

The solution of Equation (]]]) found in the first approx- 
imation by means of the Chapman-Enskog method is 
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where f\ (x a ; t) is the local quasi-equilibrium Maxwell 
one-particle distribution function: 
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This function is the solution of Equation (|l|) in the ze- 
roth approximation and satisfies the Fredgolm conditions. 
Correction ip(xf;t) reads: 
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In other words, (v) is nothing but hydrodynamical veloci- 
ty V(r;t). Functional A Q [(Cf ) 2 ], B Q [(C? ) 2 ], ^(C?) 2 ] 
are defined by the Sonine-Laguerre polynomials ||[§. 

Having the solution of Equation (|l|) in the first approx- 
imation, one can calculate the stress tensor n, heat flow 
vector q and diffusion velocity V d in the same approxi- 
mation. The expression for l~l reads: 

n = pi - k (V • V) I - 277S, 

where p is total pressure in the first approximation (I is 
the unit tensor, S is the velocities shift tensor): 

2 M 
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(g2^(x) is the binary quasi-equilibrium correlation func- 
tion, $ a a(x) is a long-range potential of interaction); k is 
the bulk viscosity of a multicomponent mixture: 
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rj is the shear viscosity of a multicomponent mixture: 
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The expression for heat flow q reads: 

M 

q=-AVT+^w tt d tt . 

a=l 

Quantities uo a are connected with a matter transfer due 
to a temperature gradient (the Soret effect) and due to a 
heat transfer caused by a gradient of concentration (the 
Dufour effect). If one also takes into account a barrodif- 
fusion process, this constitutes the total contribution into 



a heat flow from cross transfer processes. A is the thermal 
conductivity coefficient of a multicomponent mixture: 
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It should be noted, however, that not A, but the heat 
conductivity coefficient x is measured experimentally. It 
is mutually connected with A by the relation X = A/ (pC p ), 
where C p is the heat capacity at constant pressure. 
Diffusion velocity in the first approximation reads: 



where 
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is the thermal diffusion coefficient of a mixture, while 
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is the mutual diffusion coefficient. Quantities B a (0) in (||), 
A a (Q) and A a (l) in (@) and (|), E a/3 (0) in (§) are nothing 
but coefficients of expansion in the Sonine-Laguerre poly- 
nomials. Their general definition for an arbitrary potential 
of interaction is given in || . The calculations for a special 
case of a two-component mixture is performed in J/j. All 
the obtained quantities B a (0), A a (0), A a (l) and E a p(0) 
ultimately depend on the so-called 17-integrals. 

Numerical calculation for transport coefficients k, 77, 
A, and D a @ has been performed. For two- and three- 
component mixtures of neutral and charged hard spheres 
we studied the dependences of transport coefficients on 
density, temperature, and concentration ratio of some mi- 
xture components |L8| . There are a lot of approaches which 
allow to calculate viscosity with sufficient accuracy. How- 
ever, these approaches do not allow to calculate well ther- 
mal conductivity for dense and moderately dense systems. 
Our theory is devoid of such a circumstance. This is rea- 
ched due to the following: firstly, the Enskog-Landau ki- 
netic equation is obtained by means of the non-equilibrium 
statistical operator method from the first principles of 
statistical mechanics without phenomenological assump- 
tions; secondly, kinetics and hydrodynamics in the stud- 
ied systems are considered simultaneously. The last fac- 
tor is very important for dense and moderately dense 
systems (l3[|TJ,|rj,^6|,|TJ. As a result, we obtain a good 
agreement between the theory and experimental data. In 
Fig. [I] we present the concentration ratio dependences of 
some transport coefficients for two- and three-component 
mixtures for charged hard spheres. In all cases, the total 
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Fig. 1. Concentration ratio dependences of transport coefficients of some binary and ternary systems of neutral and charj 
spheres. Thermal conductivity (a) and thermal diffusion coefficient (b) of a system Ar + -Kr + at total concentration n = 
cm -3 . Thermal conductivity of a system He-Kr-Ar + (c) at total concentration n = 2 ■ 10 21 cm -3 , xat = 0.2 or tia? = 
cm -3 . Thermal conductivity of a system He-Ar + -Kr + (d) at total concentration n = 1.25 ■ 10 21 cm -3 , ihc = 0.6. 
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concentration is considered to be constant. Fig. [j] (b) il- 
lustrates the limiting cases for the thermal diffusion co- 
efficient of a system Ar + -Kr + when xkt — ► or 
XKr — ► 1) that is for one-component systems. In this case, 
diffusion thermodynamic forces vanish and thermal dif- 
fusion vanishes too. Fig. || shows numerical calculations 
for the same systems as in Fig. |[ but for point-like neu- 
tral and charged particles. Transport coefficients for last 
ones were calculated as for the usual Boltzmann-Landau 
kinetic equation ]iT| . The magnitude of transport coeffi- 
cients therewith slows down. It is interesting to note that 
in three-component systems of point-like particles (Fig. || 
(c) and (d)) transport coefficients depend on A slightly. It 
is well to bear in mind that this takes place if at least one 
component is not charged. 

From the present letter, one can draw the following 
conclusions. The obtained Enskog-Landau kinetic equa- 
tion for charged hard spheres turned out to be very useful 
for several purposes. First of all, the collision integral of 



this equation does not contain a divergency at small dis- 
tances. Secondly, the normal solution and all transport 
coefficients have analytical structure. They can be easily 
used to study some specific systems. Finally, the analyt- 
ical structure of transport coefficients allows us to find 
fast and easily systems, which can be best described by 
the obtained kinetic equation, as well as density and tem- 
perature ranges, where the agreement between the theory 
and experimental data is the closest. At the same time, our 
theory has not met with success in treatment of the free 
electrons role. Basing on the conclusions by Ichimaru et. 
al. jl9), consistent treatment of electrons is possible only 
within the frame of quantum kinetic theory. Our theory 
is purely classical. Just the same, our papers show 
very good agreement between theoretical calculations and 
experimental data. 

The next step in this theory is to calculate a dynamical 
screening radius in a system. Partially this problem has 
been already solved in our recent paper |2Cj|] . 
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Fig. 2. The same as in Fig. 1, but for systems of neutral and charged point-like particles at T = 5000 K for (a), (b) and (c) 
and at T — 6500 K for (d). Parameter A indicates for what value of Coulomb logarithm calculation was performed. 
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